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ABSTRACT 

Context. Different environmental conditions can play a crucial role in determining final products of the star formation process and in 
this context, less favorable activities of star formation are expected in the external regions of our Galaxy. 

Aims. We studied the properties of the young open cluster NGC 1893 located about 12Kpc from the galactic center, to investigate 
how different physical conditions can affect the process of star formation. 

Methods. By adopting a multiwavelength approach, we compiled a catalog extending from X-rays to NIR data to derive the cluster 
membership. In addition, optical and NIR photometric properties are used to evaluate the cluster parameters. 

Results. We find 415 diskless candidate members plus 1061 young stellar objects with a circumstellar disk or class II candidate mem- 
bers, 125 of which are also H„ emitters. Considering the diskless candidate members, we find that the cluster distance is 3.6±0.2kpc 
and the mean interstellar reddening is E(B-V)=0.6±0.1 with evidence of differential reddening in the whole surveyed region. 
Conclusions. NGC 1 893 contains a conspicuous population of pre-main sequence stars together with the well studied main sequence 
cluster population; we found a disk fraction of about 70% similar to that found in clusters of similar age in the solar neighbour and 
then, despite expected unfavorable conditions for star formation, we conclude that very rich young clusters can form also in the outer 
regions of our Galaxy. 

Key words. - - star formation 



1. Introduction 

The sequence of events that characterize the star formation 
(SF) process starts from the compression of gas in a giant 
molecular cloud that produces molecular cores; these collapse 
into protostars evolving in pre-main sequence (PMS) objects 
through accretion of mater ial and circumstellar disk formation 
dZinn ecker & YorkeL I2007I) . Within this standard picture, how- 
ever, different environmental properties may play a key role by 
changing the physical conditions during the SF process and pos- 
sibly the final products, for example the initial mass function 
(IMF). 

SF regions with peculiar properties are therefore crucial for 
a complete understanding of the SF process since they are the 
natural places where different initial conditions can lead to the 
formation of clusters with non-standard IMF, dynamical evolu- 
tion and/or disk fraction and evolution. In this context, clusters 
in the outer Galaxy are very interesting since they are located 
in regions where surface and volume densities of atomic and 
molecular hydrog e n are much smaller than in the inner Galaxy 



tion (Vall enari et al l ll999HSharma et al.Ll2007l) . This SF region, 
located at the center of the Aur OB2 associations, is associated 
with the HII region IC410 and the two pennant nebulae, Sim 19 
and Sim 130 dGaze & ShainL[l952l) . It is obscure d by several dus t 
clouds and contains at least five O-type stars dHiltnerl Q966). 
The presence of several O-type s tars and emission line B-type 
stars (M arco & Negueruelal 12002). likely in the PMS phase, is a 
strong indication that very recent SF events have occurred in this 
cluster. 

Several photometric and low resolution spectroscopic works 
have been devote d to this cluster since the first study by 
Hoag et aTl dl96ll) . Many of them are, however, aimed at 
studyin g the NGC 1893 massi ve stars, that are in the MS 
phase dBecker & FenkartL 1971 
Humphreys! 1 19781 
Mas sev et all 119951 



'l Olapia et al 
9l iMarco et al 




19721; ICuffeyL 
Fitzsimmons, 


1973 


1993 


Marco & Negueruela 



( Wou terloot et alt 1 19901). while metal content is, on average, 



smaller dwilson & Matteuccj 1 19921) and therefore higher tem 
peratures of the cloud are expected due to lower radiative losses. 
In addition, the lack of prominent spiral arms and the presence of 
few supernovae as external triggers for SF are further indications 
of less likely SF conditions. 

We have identified NGC 1893 as a very promising clus- 
ter for this kind of investigations since it is very young (~1- 
2Myr), distant about 12Kpc from the Galactic Center, with 
known massive members and indication of a large PMS popula- 



has been found by Vallenari et al.l d 19991) by using near-infrared 
photometry, even if their data did not allow them to distinguish 
the youn g population from th e contaminating field stars. Only 
recently, Shar ma et al.l (|2007j) presented a study about the low 



mass population in NGC 1893 based on deep optical photome- 
try down to V-22 and NIR 2MASS dCutri et all 120031) catalog. 
Several conclusions are drawn in this work about the IMF, the 
dynamical evolution and the effects of massive cluster stars on 
the low mass population. However, their results are based on a 
statistical subtraction of the field star population falling in the 
same region and suffer from the large statistical uncertainties re- 
lated to the membership. In fact, the main challenge of the anal- 
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ysis of the low mass population in very young clusters is the 
membership derivation since PMS stars usually lie in the same 
region occupied by contaminating field stars and it is therefore 
crucial an accurate analysis to identify cluster members. 

We have started a project aimed at the detailed study of 
the star formation process in NGC 1893 and in particular the 
unknown low mass population: the first results, based on the 
joint Chandra-Spitzer large program The Initial Mass Function 
in the Outer Galaxy: the star fo rming region NGC1893 (P.I. G. 
Micela), have been published in ICaramazza et al.l (120081 here- 
after Paper I). In this latter paper, a conspicuous sample of 
low mass cluster members with circumstellar disk and of more 
evolved diskless candidate members belonging to this region has 
been identified by using Spitzer-IRAC and Chandra X-ray data. 

We present here the results on the membership and clus- 
ter parameters based on a multiwavelength approach including 
new deep optical and JHK photometry that allow us to derive a 
more complete census of the low mass population in this cluster; 
knowledge of individual members is used to derive cluster pa- 
rameters (reddening and distance) in a more accurate way than 
previously reported in the literature. Detailed studies on the X- 
ray coronal properties, disk fraction and evolution and IMF will 
be presented in forthcoming papers. 

The present paper is organized as follow. In Sect. 2 we 
present the observations, the data reduction method and a de- 
scription of the procedure used to derive the photometry and 
the astrometry. The whole catalog is presented in Sect. 3 while 
Sect. 4 includes the comparison with previous works; Sect. 5 de- 
scribes how we identify cluster candidate members by using 
VRIJHK, Hq, and X-rays data, while in Sect. 6 the cluster pa- 
rameters are derived and a comparison with literature results is 
presented in Sect. 7. In Sect. 8 individual masses and ages are de- 
rived while in Sect. 9 we present our summary and conclusions. 

2. Observations, data reduction and photometry 

We present here new optical and near infrared (NIR) data of the 
cluster, collected from the Telescopio Nazionale Galileo (TNG) 
at the Roque de los Muchachos Observatory (ORM, Canary 
Islands, Spain), and from the 2.2 m telescope at the Calar Alto 
Observatory (Spain). 

2.1. Optical data 

Optical observations were acquired in service mode using two 
different telescopes, with the standard VRI and Ha photo- 
metric filters: the Device Optimized for the LOw RESolution 
(DOLORES) mounted on the TNG was used in service mode 
during three nights in 2007, and the Calar Alto Faint Object 
Spectrograph (CAFOS), mounted on the 2.2 m telescope in Calar 
Alto German-Spanish Observatory (Spain), observed for three 
nights in 2007 and 2008, as detailed in Table Q] Standard fields 
were taken with both instruments to perform the required photo- 
metric calibration. 

DOLORES is a focal reducer instrument equipped with a 
2048x2048 square pixels CCD, with a scale of 0.252"/px and 
a field of view (FoV) of 8.6'x8.6'. A mosaic of four fields was 
made, covering an approximate area of 22.2'xl8.0'. The stan- 
dard VRI filters were used in the four fields, as well as the Ha 
narrow band filter centered on 6658 A, with a Full Width Half 
Maximu m (FWHM) o f 57 A. We also observed the standard field 
SA98 (lLandoltL[T992l a=6 h :52 m :04.1 s , 5=-00 d :19 m :37.9 s ). 

CAFOS is a focal reducer instrument that works with a 
CCD of 2048x2048 square pixels, and has a scale of 0.53'7px, 



covering a circular FoV of 16' diameter. With this FoV we 
covered the central areas of the cluster not observed by the 
four DOLORES fields. Standard VRI filters were used, to- 
gether with an Ha filter centered at 6577 A with FWHM 
of 98 A. The fields used as standard for the photometric 
calibration were RU 149 (a=07 h :24 m :40.8 s , 5=-00 d :34 m : 10.0 s ), 
PG2331 (a=23 h :34 m :14.1 s , £=+05 d :46 m :25.9 s ), and SA95 (a= 
03 h :53 m :32.9 s , 5=+00 d :00 m :48.3 s ). In this case we noticed some 
variability in the photometric conditions of the last two nights in 
Calar Alto. To improve the quality of the photometric calibration 
we used also a set of local standards of the cluster, selected dur- 
ing the TNG data calibration, as explained below. A combined 
image of the four TNG fields in the I band, with the shape of the 
CAFOS FoV over plotted, is shown in Fig.Q] 

Optical data were reduced using the standard Image 
Reduction and Analysis Facility (IRAF) routines, including the 
bias and flatfield subtraction. Bias correction was applied by 
using Zero images and the local overscan strip in each image. 
Some non-standard steps were taken due to various problems 
arisen during in the reduction, with details described below. 

For the DOLORES observations a combination of sky flat- 
fields was used to create a master flat field in each band, using 
the task flatcombine. Once the images were corrected from flat- 
field, a strong illumination pattern was still present in the images 
of the I band. An illumination function was created using the flat- 
field (with the task mkillumflat), and conveniently removed from 
images. A pattern of interference fringing was still present in the 
I band images, with no variations found in the different nights. 
We constructed a fringing image by combining 15 images in the 
I band of the three nights, so the stars in the different fields are 
canceled. The images selected were of short exposure (2-60 s) 
to reduce the number of stars in the fields. The resulting fringing 
pattern was successfully corrected using the task mkfringecor 
and applied to every image in the I band. 

For the CAFOS data, dome flats were used to correct the im- 
ages from local irregularities. The singular circular shape of the 
FoV required special care: the values of all the pixels outside the 
FoV (easy to filter because they are well below the background) 
were set to a very large number; in this way no unreal features 
are created in the target images after subsequent flat-field re- 
moval, resulting in a value well below background outside the 
FoV. The remaining reduction followed the standard path. 

2.2. NIR data 

NIR observations were acquired in service mode at the TNG, us- 
ing the large field Near Infrared Camera Spectrometer (NICS), 
that is a camera equipped with a 1024x1024 square pixels 
HgCdTe HAWAI I infrared array, w ith a scale of 0.25"/pix and a 
FoV of 4.2'x4.2' dBaffaet all 1200 lb . Our observations were ac- 
quired with the Js(1.25ju), H(1.63//) and K'(2.12ju) filters dur- 
ing eight nights in 2007 and 2008; we observed the whole field 
with a raster of 4x4 pointings, at each pointing we obtained a 
series of NINT dithered exposures; each exposure is a repetition 
of a DIT (Detector Integration Time) times NDIT (number of 
DIT), to avoid saturation of the background; details are given in 
Table|2 

Raw IR images taken with NICS were reduced by using the 
Speedy Near-IR data Automatic reduction Pipeline (SNAP) im- 
plemented for NICS using the "full dither" options. We forced 
the script to use an external master flat field obtained by a median 
combination of single flat fields taken during the same night or 
the nearest night. After a bad pixel correction, the pipeline per- 
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Fig. 1. Upper panel: I band combined image of the four DOLORES frames of 740 sec; the overplotted circle indicates the CAFOS 
FoV. Bottom panel: Js band combined image of the NICS frames reduced and coadded with SNAP. 
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Table 1. Optical observations of the cluster NGC 1893 



Date 



a (J2000) 5 (J2000) Instrument Filter Exp. time[sec] seeing airmass 



2007/09/21 5:22:27.4 33:23:46.0 DOLORES V 10,60,2x500 1.0" 1.151-1.192 

2007/09/21 5:22:27.5 33:23:36.9 DOLORES R 10,70,700 1.0-1.1" 1.066-1.079 

2007/09/21 5:22:27.4 33:23:42.5 DOLORES I 10,60,2x740 0.9" 1.095-1.138 

2007/09/21 5:22:27.6 33:23:33.0 DOLORES Ha 60,300,2x700 0.9-1.0" 1.024-1.049 

2007/10/18 5:23:11.6 33:23:34.3 DOLORES V 10,60,2x 500 0.8" 1.026-1.040 

2007/10/18 5:23:11.9 33:23:37.1 DOLORES R 10,70,700 0.8-1.0" 1.006-1.010 

2007/10/18 5:23:11.8 33:23:22.0 DOLORES I 10,60,2x740 0.8-1.0" 1.009-1.021 

2007/10/18 5:23:12.8 33:23:47.2 DOLORES Ha 60,300,2x700 0.8" 1.013-1.031 

2007/11/14 5:23:06.9 33:32:29.7 DOLORES V 10,60,2x 500 1.0" 1.282-1.353 

2007/11/14 5:23:06.8 33:32:14.6 DOLORES R 10,70,700 1.0-1.1" 1.113-1.134 

2007/11/14 5:23:06.7 33:32:24.5 DOLORES I 10,60,2x740 0.8-0.9" 1.184-1.257 

2007/11/14 5:23:06.7 33:32:08.2 DOLORES Ha 20,60,300,2x700 0.8" 1.050-1.096 

2007/11/14 5:22:28.7 33:32:30.5 DOLORES V 10,60,2x500 1.0" 1.004-1.007 

2007/11/14 5:22:33.9 33:32:39.3 DOLORES R 2x 10,70,700 0.8-0.9" 1.015-1.022 

2007/11/14 5:22:32.0 33:32:20.5 DOLORES I 10,60,2x740 0.9-1.0" 1.004-1.010 

2007/11/14 5:22:34.2 33:32:47.4 DOLORES Ha 60,2x 300,2x700 0.8" 1.029-1.066 

2007/10/11 5:23:22.0 33:26:38.7 CAFOS V 15,3x 500 1.4-1.6" 1.002-1.004 

2007/10/11 5:23:21.5 33:26:37.4 CAFOS Ha 15,100,3x 500 1.6-1.8" 1.002-1.016 

2008/01/05 5:23:23.3 33:26:42.4 CAFOS V 10,60,2x300 2.4" 1.161-1.205 

2008/01/05 5:23:23.2 33:26:45.8 CAFOS R 10,60,2x 300 2.3-2.6" 1.231-1.283 

2008/01/05 5:23:23.3 33:26:47.3 CAFOS I 10,60,3x 500 1.5-1.6" 1.353-1.502 

2008/01/05 5:23:23.1 33:26:38.9 CAFOS Ha 10,60,500 1.4" 1.603-1.651 

2008/01/09 5:23:24.2 33:26:41.6 CAFOS Ha 15,100,3x 500 1.2-1.7" 1.139-1.239 



forms the source detection on the sky subtracted and flat fielded 
images. A double iteration is performed to obtain a coadded im- 
age for each filter by using the offsets among the images com- 
puted by a cross-correlation algorithm after the presence of the 
field distortion is corrected. A combined image in the Js band of 
the NICS images reduced and coadded with SNAP is shown in 
Fig. Q] (bottom panel). 



2.3. Photometry and data selection 

Sources detection and Point Spread Function (PSF) photome- 
try were made using the routines within the Fortran sta nd-alone 
version s of the DAOPHOT II and ALLSTAR pac kages (IStetsonL 
Il987ah . improved with the use of ALLFRAME (Stetson Ul99ll) . 
Instrumental PSF photometry was corrected by using appropri- 
ate aperture correction s derived by growth curves computed with 
the DAOGROW code dStetsoij[l990l) . Details on the photomet- 
ric calibration are given in Appendix. 

To perform an accurate analysis of the color-magnitude and 
color-color diagrams and compare them with appropriate the- 
oretical models we selected the optical photometric catalogs 
obtained from the Dolores and Cafos images by considering 
only objects with photometric error smaller than a magnitude- 
dependent limit and with the DAOPHOT parameter sharp 
smaller than 0.5 i.e. objects with brightness distribution consis- 
tent with point-like sources. The Dolores and Cafos final cata- 
logs include 7144 and 3222 objects, respectively. 

The photometric JHK catalog obtained using NICS images 
was selected by considering only objects with magnitude errors 
smaller than ^4 in order to have color errors smaller than 0.2. 

V2 

This criterion was adopted in order to discard false detections 
found in the spots due to saturated objects. After this selec- 
tion our JHK catalog contained 11181 entries including multiple 
items for the objects located in the overlapping regions of the 16 
fields. 



Offset [arcsec] 



Fig. 2. Cumulative distribution of offsets between objects falling 
in adjacent fields (dotted line). The expected offset cumulative 
distribution of spurious matches (dashed line, see Eq.[TJ and the 
true distribution (solid line, see text) are also shown. 



To filter out the redundant detections, we first estimated an 
appropriate matching radius by computing the offset cumulative 
distribution of the objects in adjacent fields within a relatively 
large radius equal to 3" using a binsize of 0.2", as shown in 
Fig.|2](dotted line). 

We compared this distribution with the expected cumula- 
tive distribution of spurious identifications (dashed line in Fig.|2j 
computed by using the relation 



7»r 15,16 

dN _ y> ii, * // 

17 ~ 2-i 



(=1,7=2 Ai J 



i-[n{r + drf - nr 2 } 



(1) 



where n,- and rij are the total number of objects in the i-th 
and j-th fields, respectively, falling in the common area Ay. The 
difference between the measured offset distribution and the ex- 
pected spurious one is the true distribution (solid line in Fig.|2j 
which tends to flatten at separations greater than 1.1". 
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Table 2. Near Infrared observations of the cluster NGC 1893 from NICS; NINTs is the number of single frames. 
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Date 


Field 


a 


(J2000) 


6 (J2000) 


Filter Exp. time [sec] NINTs 


seeing 


airmass 


2007/10/08 


10 


5: 


:22:59.1 


33:25:20.2 


H 


600 


15 


0.8" 


1.103 


2007/10/08 


10 


5 


:22:59.0 


33:25:23.3 


J 


510 


17 


0.8" 


1.133 


2007/10/08 


10 


5: 


:22:58.8 


33:25:26.6 


K' 


700 


14 


0.8" 


1.178 


2007/10/08 


5 


5: 


:22:46.9 


33:21:48.1 


H 


600 


15 


1.0" 


1.003 


2007/10/08 


5 


5 


:22:44.4 


33:21:10.7 


J 


510 


17 


1.2" 


1.004 


2007/10/08 


5 


5 


:22:41.5 


33:21: 1.1 


K' 


700 


14 


0.8" 


1.009 


2007/10/08 


6 


5: 


:22:42.6 


33:25:15.6 


H 


600 


15 


1.0" 


1.022 


2007/10/08 


6 


5: 


:22:42.1 


33:25:20.0 


J 


510 


17 


1.0" 


1.035 


2007/10/08 


6 


5 


:22:41.6 


33:25:24.7 


K' 


700 


14 


1.0" 


1.056 


2007/10/10 


11 


5 


:22:59.3 


33:29:15.7 


H 


600 


15 


0.8" 


1.197 


2007/10/10 


11 


5: 


:22:59.1 


33:29:18.6 


J 


510 


17 


1.0" 


1.243 


2007/10/10 


11 


5: 


:22:59.0 


33:29:22.6 


K' 


700 


14 


0.8" 


1.314 


2007/10/10 


9 


5:23: 0.2 


33:21: 2.8 


H 


600 


15 


1.0" 


1.073 


2007/10/10 


9 


5 


:22:60.0 


33:21: 6.0 


J 


510 


17 


1.0" 


1.096 


2007/10/10 


9 


5: 


:22:59.8 


33:21: 9.8 


K' 


700 


14 


1.0" 


1.133 


2007/10/19 


12 


5: 


:22:59.6 


33:33:59.9 


H 


610 


16 


0.8" 


1.467 


2007/10/19 


12 


5 


:22:59.7 


33:34: 0.0 


J 


510 


17 


0.8" 


1.375 


2007/10/19 


12 


5 


:22:59.7 


33:34: 0.0 


K' 


700 


14 


0.8" 


1.313 


2007/10/19 


15 


5: 


:23:17.2 


33:28:55.8 


H 


600 


15 


0.9" 


1.186 


2007/10/19 


15 


5: 


23:17.3 


33:28:53.7 


J 


510 


17 


1.0" 


1.144 


2007/10/19 


15 


5 


:23:17.5 


33:28:50.5 


K' 


700 


14 


1.0" 


1.114 


2007/10/19 


16 


5 


:23:19.0 


33:32:58.1 


H 


600 


15 


1.0" 


1.034 


2007/10/19 


16 


5: 


:23:19.0 


33:32:58.1 


J 


510 


17 


0.9" 


1.049 


2007/10/19 


16 


5: 


:23:18.7 


33:34: 5.0 


K' 


710 


15 


0.9" 


1.078 


2007/10/19 


7 


5 


:22:40.7 


33:30: 5.1 


H 


600 


15 


0.9" 


1.571 


2007/10/19 


7 


5 


:22:40.6 


33:30: 5.1 


J 


510 


17 


0.9" 


1.838 


2007/10/19 


7 


5: 


:22:40.7 


33:30: 4.8 


K' 


700 


14 


0.9" 


1.707 


2007/12/17 


1 


5: 


:22:21.9 


33:21:15.0 


H 


600 


15 


1.2" 


1.992 


2007/12/17 


1 


5: 


:22:21.9 


33:21:15.0 


K' 


700 


14 


1.2" 


2.272 


2008/01/01 


3 


5 


■22:21.1 


33:30:52.3 


H 


600 


15 


1.0" 


2.047 


2008/01/01 


3 


5 


:22:27.6 


33:30:52.5 


J 


510 


17 


1.0" 


1.882 


2008/01/01 


3 


5: 


:22:21.1 


33:30:52.3 


K' 


700 


14 


1.2" 


2.290 


2008/01/11 


1 


5: 


:22:23.9 


33:21:53.5 


J 


510 


17 


1.2" 


1.419 


2008/01/11 


4 


5 


:22:23.9 


33:33:53.5 


J 


510 


17 


1.5" 


1.642 


2008/01/13 


4 


5 


:22:30.2 


33:34: 5.0 


H 


650 


16 


1.0" 


1.019 


2008/01/13 


4 


5: 


:22:30.2 


33:34: 4.9 


fC 


700 


14 


1.0" 


1.034 


2008/01/16 


13 


5: 


:23:25.9 


33:22:22.1 


H 


600 


15 


0.8" 


1.151 


2008/01/16 


13 


5 


:23:25.9 


33:22:22.1 


J 


510 


17 


0.8" 


1.196 


2008/01/16 


13 


5 


:23:25.9 


33:22:22.1 


K' 


700 


14 


0.9" 


1.248 


2008/01/16 


14 


5: 


:23:24.8 


33:26:43.0 


H 


600 


15 


1.0" 


1.586 


2008/01/16 


14 


5: 


:23:24.8 


33:26:42.9 


J 


510 


17 


0.8" 


1.489 


2008/01/16 


14 


5 


23:18.6 


33:26: 5.8 


K' 


760 


16 


0.8" 


1.342 


2008/01/16 


2 


5 


:22:20.4 


33:25:10.2 


H 


138 


31 


1.2" 


1.170 


2008/01/16 


2 


5: 


:22:20.4 


33:25:10.2 


J 


540 


18 


1.2" 


1.213 


2008/01/16 


2 


5: 


:22:18.4 


33:25:38.5 


K' 


700 


14 


1.0" 


1.047 


2008/01/16 


8 


5 


:22:39.8 


33:32:47.8 


H 


560 


14 


1.2" 


1.006 


2008/01/16 


8 


5 


:22:39.8 


33:32:47.8 


J 


510 


17 


1.0" 


1.011 


2008/01/16 


8 


5: 


:22:39.8 


33:32:47.8 


K' 


700 


14 


1.0" 


1.023 



Using this value as the final matching radius, we included 
less than 10% of spurious matches and most of the true identifi- 
cations. 

Therefore, we matched each star of our catalog with those of 
adjacent fields and considered as the same object, those falling 
within 1.1". The final catalog composed of single items includes 
9 880 objects with K magnitude down to 17.78. 

2.4. Astrometry 

In order to assign celestial coordinates to the stars in the opti- 
cal catalogs, we used the 2MASS catalog dCutri et all 120031) as 
reference. The transformation between the two coordinate sys- 



tems was made by applying the appropriate sky projection us- 
ing several tasks of the imcoords IRAF package. To determine 
this transformation we first identified through visual inspection 
three stars common to the 2MASS observations and to each of 
our fields. With these initial coordinates we used the task ccxy- 
match to match the pixel coordinate list of each of our fields 
and the celestial coordinate list of the 2MASS catalog. The re- 
sulting single matched coordinate list is used as input of the task 
ccmap to derive the final sky projection. This sky map is set 
in the master image of each field with the task CCSetwcs, and 
the catalog of stars converted into celestial coordinates using the 
xy2sky program of the WCSTools UNIX package developed 
at the Smithsonian Astrophysical Observatory. The assignment 
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of celestial coordinates to the stars in the NIR catalog was per- 
formed by using the WCSTools package. 

To verify the results and estimate the astrometric accuracy, 
we matched the three entire catalogs DOLORES, CAFOS and 
NICS, with the 2MASS catalog, used as reference to find the as- 
trometric solution, and considered the offset distribution within 
a relatively large value (2"). From this distribution we sub- 
tracted the expected distribution of spurious identifications (Eq. 
1), and we obtained the distribution of true identifications. For 
the matches DOLORES -2M ASS, CAFOS-2MASS and NICS- 
2MASS, the 1 <x corresponding values are 0.24", 0.22" and 
0.19", respectively, that are then the final accuracies for each 
catalog. 

3. Catalog 

The aim of this section is to describe the multiwavelength cata- 
log obtained using the following data: 

- the new optical VRI photometry from DOLORES® TNG 
and CAFOS @2. 2m telescope of the Calar Alto Observatory 
presented here; 

- the new deep JHK photometry from NICS @ TNG telescope 
presented here; 

- the Spitzer/IRAC photometry published in Paper I; 

- the Chandra/ ACIS -I X-ray detections published in Paper I. 

Figure[3]shows the fields of view (FoV) of all involved obser- 
vations, viz. DOLORES (4 solid boxes), CAFOS (solid circle), 
NICS (16 dotted boxes) plus the Chandra-ACISQ (dashed box). 
To match these catalogs we first derived an appropriate match- 
ing radius for each couple of catalogs, and then we identified 
the common objects by considering the nearest neighbor match 
within the adopted radius. With the aim to include as many coun- 
terparts as possible between the catalogs, we initially computed 
the offset distribution within a radius larger than the expected 
astrometric accuracy of the two matched catalogs and the dis- 
tribution of the spurious identifications according to the Eq. Q] 
To match the objects we use the radius at which the true distri- 
bution (see Sect. 12.3b starts to flatten; in all cases we have less 
than 10% of spurious identifications. We derived a matching ra- 
dius of 0.7" between the NICS and IRAC catalogs. The IRAC 
catalog includes 15 277 objects having at least one magnitude 
among the four IRAC filters; 7194 sources fall within the NICS 
FoV, 5020 with a counterpart in the NICS catalog. Most of the 
NICS objects with no IRAC counterpart are fainter than K=15, 
since NICS data are deeper than IRAC observations. 

Using the same method described before, we found that 
a matching radius of 0.7 " is appropriate to cross-match the 
DOLORES optical catalog with both the NICS and IRAC 
catalogs. In the case of the DOLORES-NICS matches we 
found 5285 common objects; using the same radius, we cross- 
correlated the DOLORES objects without a NICS counterpart 
with the IRAC catalog and we found further 629 objects, mostly 
located outside the NICS images, for a total of 5914 DOLORES 
optical detections having at least one NICS or IRAC counterpart. 
A radius of 0.7 " was used also to cross-correlate the CAFOS 
and NICS catalogs, finding 2915 common objects, while a ra- 
dius of 0.9" was adopted to match the CAFOS unmatched ob- 
jects with the IRAC catalog, finding further 64 objects. 

The X-ray ACIS-I catalog includes 1021 sources. As men- 
tioned in Paper I, 717 of them have a single counterpart in the 
IRAC catalog while 2 sources have a double identification. For 

1 Note that the IRAC FoV is larger than the fields shown in this figure. 



Table 3. Comparison between our and literature optical and NIR 
photometric catalogs. For each couple of catalogs, statistical val- 
ues and the number (N) of objects on which they have been com- 
puted are given. 
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the match between the NICS and X-ray ACIS-I catalogs we used 
the same procedure adopted in Paper I, i.e. we considered three 
different matching radii, 0.7", 1.5" and 2.5", in order to take 
into account the degrading of the ACIS point spread function 
(PSF) off axis. We found a total of 752 X-ray sources with a sin- 
gle counterpart in the NICS catalog and 10 X-ray sources with 
two IR counterparts 642 having already an IRAC detection. 

The entire catalog includes 21 688 objects, with 9531 having 
only the IRAC counterpart since the Spitzer/IRAC FoV is larger 
than that of the other observations. 



4. Photometry comparison 

We compared the photometry obtained with DOLORES and 
CAFOS data in order to check the consistency of our data. 
We limited the comparison to V=22 since the limiting magni- 
tudes in the CAFOS and DOLORES catalogs are V~22 and 
V~ 24.3, respectively, if we consider magnitudes with errors 
smaller than 0.1. The results of the statistical comparison be- 
tween the catalogs are shown in Table[3]where we also show the 
an alogous result s obtain ed from the comparison with the catalog 
of ISharma et alJ (120071) . that is the only VRI photometric cata- 
log avaikdbJe_mHterature. For the latter comparison, we limited 
the ISharma et al. I J2007I) catalog to V<20 in order to avoid ef- 
fects due to large photometric errors. Our catalogs, obtained with 
CAFOS and DOLORES observations are in excellent agree- 
ment, as expected, since the CAFOS one has been calibrated by 
using the DOLORES c atalog . We found a good agreement also 
with the Sharma et al. (2007) photometry since in all bands we 
have offsets smaller than 0.04 mag and the standard deviations is 
always smaller than 0.2 mag. The mean offsets in the V-I color 
are all smaller than 0.05 mag. 

A good agreement is found also from the comparison be- 
tween our JHK photometry, obtained with NICS data, and the 
2MASS catalog, as shown by the corresponding statistical val- 
ues given in Table[3] 

5. Member selection 

In this section we complement our optical phot ometric cata- 
log by using the literature V and I magnitudes bv lSharma et al.l 
(2007) for the objects detected in X-rays with Chandra/ ACIS -I 
or in the NIR with NICS and/or IRAC falling outside our optical 
fields or in saturated or corrupted regions of the V and I images. 
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Fig. 3. FoVs of our observations: DOLORES (solid boxes), CAFOS (solid circle), NICS (dotted boxes) plus the Chandra- ACIS 
(dashed box). In the upper-left and bottom-left corners of the DOLORES (D) and NICS (N) boxes, respectively, the corresponding 
field number is indicated. The IRAC FoV is larger than all the previous ones and it is not shown here. 



This is only for the purpose to find an as complete as possible list 
of cluster candidate members, while in the next sections, where 
we derive cluster parameters, stellar masses and ages, we will 
use only our photometry to be sure to have an homogeneous 
photometric system. 



5. 1. Candidate members with circumstellar disk 



published in Paper I, based on X-ray detection and/or infrared 
photometry in the four IRAC bands. 

As in iDamiani et all (120061) and in iGuarcello et all (120071 

2009) we selected YSOs with circumstellar disk by using 
reddening-independent indices that allow us to distinguish ob- 
jects reddened by the interstellar medium from those with IR 
excesses due to the disk. We have defined a generic index as 



Deep infrared observations in the JHK bands presented in this 
paper allow us to enlarge the sample of NGC 1893 members 



Gabcd = (A - B) 



E(A - B) 
E(C - D) 



x (C - D) 



(2) 
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Table 4. Couples of colors used to define the reddening inde- 
pendent Q-indices as in Eq.|2] and the limit adopted to distin- 
guish objects with normal photospheric colors from YSOs with 
excesses. 



A-B C-D Q [im A-B C-D Qj, 



where (A-B) and (C-D) are two different colors while E(A- 
B)/E(C-D) depends only on the adopt ed reddening law. In our 
case we used the redd e ning l aws of iMunari & Carrara ([1996) 
and [Rieke & Lebofskvl dl985h . for optical and JHK bands, re- 
spectively, since these are those that better follow our data in the 
color-color diagrams; for the IRAC bands we adopted the mean 
of the extinction s relati ve to the K band, derived for five SFRs 
in Flahert y et alj (|2P07). With the magnitudes available in our 
multiband catalog we have defined 20 reddening independent 
Q-indices by using the couples of colors given in Table|4] These 
indices allow us to select in a way as complete as possible, all 
the objects showing excesses in a given color. For example, in- 
dices involving the (V-I) colors should select YSOs with excess 
in the NIR band, if the (V-I) colors derive from the photosphere. 
However, we note that in case of accretion, the V band can be 
affected by the veiling contribution; in this case, the Q index is 
the combination of excesses in two colors. 

For each reddening independent Q-indices we adopted a con- 
servative limit that c hosen by considering the exp ected colors of 
main-sequence stars dKenvon & Hartmanni 1 19951) for the colors 
involving the VIJHK band. For the color involving the IRAC 
magnitudes [3.6], [4.5], [5.8] and [8.0] we defined the limit by 
considering the Q-indices of the objects that in the [3.6]-[4.5] vs. 
[5.8]-[8.0] diagram have color compatible with zero within 2 cr. 
By assuming these limits, given for each index in Table|4] we 
consider objects with infrared excess those with Q-index smaller 
within 3 <x than the adopted limit in a given index. 

We note however that with the Q-index method all the se- 
lected objects have a disk but there is an ambiguity region where 
stars with disk cannot be distinguished by reddened objects and 
then the list of candidate members with a disk can be incomplete. 

With this criterion, we selected 984 class II YSOs; by com- 
paring the class II YSOs with the 242 class II YSOs selected in 
Paper I we found 187 are common to the two samples, only 19 
are class II according to the IRAC colors but not according to 
our indices, and 36 are class II according to the IRAC colors that 
fall outside of the NICS Fo V. 

In our list of class II candidate members we found 5 of the 
7 YSOs classified in Paper I as class 0/1 candidate members, for 
which we maintained the 0/1 classification and discarded them 
from our sample of class II candidate members. In summary, 
we have 1034 (984+19+36-5) class II YSOs of which 792 are 
new candidate members with a disk found in this work. The lat- 
ter include 17 objects classified in Paper I as diskless candidate 
members: we checked that the NICS photometry is of good qual- 
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Fig. 4. V vs. V-I diagram of all objects with optical magni- 
tudes in our catalog. Class II YSOs are indicated wit h squares- 
Solid line is the lOMyr solar metallicity isochrone of lSiess et alj 
(2000) drawn assuming a distance of 3.6 kpc and E(B-V)=0.6. 



ity and therefore we classified these objects as class II candidate 
members. 

The V vs. V-I diagram of the objects with circumstellar disk 
is presented in Fig. [4] where all class II candidate members are 
indicated by squ ares. The solid line is the 10 Myr solar metallic- 
ity isochrone of ISiess et al.1 (|2000), drawn assuming a distance 
of 3.6 kpc and E(B-V)=0.6 (see Sect.[6]); fro m now on the the- 
oretica l m odels will be plotted by u sing the IMunari & Carrarol 
(Il996h and lRieke & Lebofskvl (1198 5) redde ning laws for th e op- 
tical and NIR bands and the iKenvon & Hartmanni (1 1995b con- 
versions to transform theoretical temperatures and luminosities 
in the observational plane. The 10 Myr isochrone is used here 
only to distinguish the PMS region from the stars that have col- 
ors apparently inconsistent with a PMS nature. 

We found that among the 1034 class II candidate members, 
there are 170 YSOs in the region below the 10 Myr isochrone. 
We do not believe these to be actually older than the other mem- 
bers, but rather, we guess they are candidate members for which 
optical colors are not purely photospheric. Most of these ob- 
jects show excesses in the JHK-IRAC colors and therefore the 
anomalous optical colors and/or magnitudes could be due to phe- 
nomena related to the disk presence, such as veiling, scattering 
or other effects due to the angle of inclination of the disk (see 
iGuarcello et aIll20Tol for a discussion). 

In fact, using the synthetic photometry derived from the 
Robi taille et al.l (120061) models, stars with circumstellar disk with 
age younger than 10 Myr can have optical colors and magni- 
tudes apparently compatible with older age. However we cannot 
exclude that a small fraction of these objects could also be ex- 
tragalactic sources, as expected since the cluster NGC 1893 lies 
in the direction of the galactic anti-center. We defer to a future 
work an accurate analysis of the disk properties of these objects 
by means of spectroscopic data. 
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Fig. 5. V vs. V-I diagram of all objects with optical magnitudes 
in our catalog. X-ray detections are indicated with x symbols 
while those selected as diskless YSOs are indicated with filled 
circles. Solid line is as in Fig. [4] 



5.3. H a emitters 

In the case of stars with accretion from circumstellar disk, the 
H a line is much enlarged and therefore the (R-H t? ) index is a 
measure of the H a flux in excess. In order to identify stars with 
H a emission, we used the R vs. R-H ff diagram shown in Fig. [6] 
where we defined the normal star limit by performing a linear fit 
of the R magnitudes of all diskless YSOs (indicated with filled 
circles) having R-H a < -3. Then we estimated a typical error 
in the R-H„ index per bin of 0.5 magnitudes, as the average of 
the R-H ff errors in each bin. We added such errors to the R-H a 
derived from the linear fit to obtain the limit indicated in the 
figure by the solid line. We defined as H a emitters, those with R- 
H a larger within 3 cr than this limit. Among the 269 objects that 
satisfy this condition we consider as cluster H a emitters only the 
125 stars (indicated in the figure as empty circles) for which we 
have an additional membership indication, viz. 98 class II and 27 
objects previously classified as diskless candidate members. 

We included in the class II sample these 27 H a emitters, that 
we had classified as diskless in the previous section. We note that 
this is not a contradiction in the classification but only a limit in 
the selection of candidate members with disk since there is a 
region in which the excesses by disks cannot be distinguished 
from reddened objects. 

In summary, we have a sample of 1061 class II candidate 
members, 125 of which are H a emitters, plus 415 diskless stars. 
Optical photometry and classification of class II and diskless 
candidates members are given in Tables|6] and [7] while new 
JHK photometry and Spitzer-IRAC magnitudes from Paper I of 
class II and diskless candidates members are given in Tables|8] 
andlH 



5.2. Candidate diskless members 

We define as diskless candidate members all the objects having 
an X-ray counterpart in the Chandra- ACIS catalog, showing any 
IR excess, i.e. not belonging to the sample of class II YSOs, de- 
fined in the previous section, and falling in the PMS region of the 
V vs. V-I diagram, i.e. that with colors redder than the lOMyr 
solar metallicity isochrone of Sies s et al ] (l2000h drawn assuming 
a distance of 3.6 kpc and E(B-V)=0.6 (see Sect.|6]l. 

To define the PMS region, we have been guided by the bulk 
of objects with X-ray counterpart, as shown in Fig. [5] which de- 
picts the V vs. V-I diagram of all objects with optical magnitudes 
in our catalog. Stars with X-ray counterpart are indicated with x 
symbols. 

In the sample of diskless candidate members we also include 
the objects with X-ray emission that do not show any IR excess 
and have already reached the MS phase at the cluster age, i.e. 
those with V<15. 

With the above conditions, we found 434 diskless YSOs, 85 
of which are in the sample of 1 10 diskless YSOs found in Paper 
I. If we discard from the 110 diskless YSOs defined in Paper 
I the 17 objects considered here as class II YSOs, the remain- 
ing 8 diskless YSOs classified in Paper I do not have either V 
or I magnitudes. By using the literature B and V magnitudes by 
iMassev et al.l d 1995b of these 8 objects, we found that 5 fall in 
the cluster MS region while the other 3 lie in the PMS region: 
we therefore included them in our sample of diskless cluster can- 
didate members. 

In summary, we have a total of 442 diskless YSOs. These 
objects are indicated with filled circles in the V vs. V-I diagram 
presented in Fig. [5] 



6. Cluster parameters 

Metallicity, interstellar reddening and cluster distance are fun- 
damental to compare the observed CMD with theoretical tracks 
and isochrones from which we estimate masses and ages of clus- 
ter candidate m embers. However, as already discussed i n sev- 
eral papers (e.g. lHillenbrandl[l997l:lGullbring et a l.. 1998), class 
II YSOs typically experience accretion phenomena, scattering 
and/or reflection processes due to the dust grains in their circum- 
stellar disks that may alter their photospheric colors, also in the 
optical bands. In such cases comparison of observed magnitudes 
and colors with models can be unreliable and the results should 
be interpreted with caution. In order to get around this problem, 
we derive cluster parameters by considering only cluster can- 
didate members without circumstellar disk (class III YSOs) for 
which the observed colors can be directly compared with evolu- 
tionary tracks and isochrones once the interstellar reddening is 
taken into account. 

6.1. Metallicity 

In this work we adopt models computed for solar metallicity, 
since there are ambiguous indications in the literature about the 
metallicity in NGC 1893. I n a pa per f ocused on the Galactic 
metallicity, IRolleston etall d 19931) and iRolleston etail {2000) 
derived C, N, O, Mg, Al and Si abundances of 80 B-type main 
sequence stars, 8 of which are NGC 1893 members. After dis- 
carding two high rotator objects, they performed an LTE analy- 

2 Full Tables[6] [7] [8] and [9] are also available in electronic 
form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr or via 
http://cdsweb.u-strasbg.fr 
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Fig. 6. R vs. R-H a diagram of all objects; diskless YSOs are in- 
dicated with filled circles while H„ emitters with at least an other 
membership indication are indicated as empty circles. Solid line 
is the limit we defined to distinguish diskless YSOs from high 
Hq, emitters. 



sis from which they derive for NGC 1893 onl y a slight indica- 
tion o f under solar abundance. More recently, Daflo n & Cunhal 
(2004) derived non-LTE chemical abundances in young OB 
stars from high-resolution echelle spectra. The two NGC 1893 
cluster members show again a marginal indication of subso- 
lar metallicity, but from the chemical analysis of 69 young 
stars with Galac t ocentr ic distances between 4.7 and 13.2kpc, 
Daflo n & Cunhal (|2004) find metallicity slopes consistent with 
a flattening of the radial gradient. Even if Rollest on et alj (2000) 
found indication of steeper slopes, the available studies do not 
show strong evidence of sub-solar metallicity and therefore, as 
assumed in previous studies, we adopt a solar metallicity for 
NGC 1893. 

6.2. Interstellar reddening 

Color-color diagram using UBV photometry is the classical tool 
to derive the mean interstellar reddening of a given cluster by 
comparing expected photospheric colors with those observed for 
diskless cluster members. Fig.Qshows the U-B vs. B-V diagram 
obtained by us ing literature UB V photometry by iMassev et al.l 
d 1995b and lSharma et al.l (120071) ; diskless stars are indicated with 
filled circles. We find that most of the bluest diskless YSOs, 
indicated by empty squares, are compatible with a reddening 
E(B-V)=0 . 6+0.1 from the c omparison with the solar metallicity 
ZAMS of Siess et al] d2000l). transformed in to the observed col- 
ors by using the lKenvon & Hartmannl(ll995l) relations. Diskless 
candidate members with B-V>0.6 can be compatible with lower 
reddenings, down to E(B-V)=0.2, or with higher reddenings, at 
least up to E(B-V)=1.3, as shown by the ZAMS reddened with 
E(B-V)=0.2 and E(B-V)=1.3 (dashed lines in Fig.|7). 

This allows us to deduce that the cluster region is affected 
by a significant differential reddening and that for the reddest 
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Fig. 7. U- B vs. B-V diagram obtained by using Masse vet al.l 
( 1995) and lSharmae t al. (2007) photometry (dots); filled circles 
are diskless candidate members, squares are those compatible 
with reddening 0.5 < E(B - V) < 0.7 based on the expected 
photospheric colors bv lKenvon & Hartmannld 19951) drawn at the 
E(B-V) values i ndicated on each line. The reddening vector ob- 
tained using 

fhe|M 

unari & Carrarol (1 19961) reddening law is also 

shown. 



objects the interstellar reddening cannot be photometrically de- 
rived. On the contrary, it can be unequivocally derived for the 
bluest objects, i.e. those indicated with squares for which the 
0.5 <E(B-V)< 0.7. This result has been checked by comput- 
ing individual reddenings for the 12 diskless YSOs with U- 
B< 0.5 for which spectral types are available dHiltnerl 1 1966; 
IMassev et aH 1 19951 Marco & Negueruelall2002l) : we found that 
among the bluest objects (B-V<0.5), all the reddening values are 
compatible with E(B-V)=0.6±0.1, while among the reddest ones 
(B-V>0.5), four have E(B-V)~0.2 and one has E(B-V)=0.68. 

The mean reddening E(B-V)=0.6±0.1 we derive for 
NGC 1893 is consistent also in the V vs. V-I CMD shown in 
Fig. [8] In this figure, dots indicate all objects in the FoV with er- 
ror in V-I smaller than 0. 1 and filled circles indicate diskless can- 
didate cluster members. Objects indicated also with squares are 
those with E(B-V) between 0.5 and 0.7, as in F ig. [7] The thick 
line is the solar metallicity ZAMS of ISiess et alj (120001) at a dis- 
tance of 3.6 kpc (see Sect. 16.31 ) and reddened using E(B-V)=0.4. 
The value E(B-V)=0.4 is that assumed for foreground objects 
that we expect to be less reddened than cl uster members. To 
fit the highest mass range not covered by the lSiess et alJ ([2000) 
ZAMS, w e use the theo r etica l isochrone at solar metallicity of 
1.5 Myr of iMarigo et all d2008l) reddened with E(B-V)=0.5, 0.6 
and 0.7 at a distance of 3.6 kpc (dashed lines). In this diagram, 
diskless candidate cluster members show a large spread that can 
be due to several reasons, such as differential reddening, binarity 
and/or age spread. This is true also in the V< 16 range where 
we can have objects in PMS or MS phase. However, those with 
0.5<E(B-V)<0.7, as derived from Fig.|7](indicated with squares) 
are well fitted by the 1.5 Myr isochrone if a reddening E(B- 
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V)=0.6±0.1 is adopted, in agreement with the value we derived 
by using the U-B vs. B-V diagram. On the contrary, the other 
bright objects show a significant spread and can be both PMS ob- 
jects or very reddened MS members. We note also that the choice 
of the 1 .5 Myr isochrone is not critical for this result since in this 
mass range, isochrones degenerate in a vertical shape branch. 
For this reason in the V vs. V-I diagram, MS members cannot 
be used to constrain neither ages, nor the cluster distance, while 
they can constrain the cluster reddening. 

Our result is also consistent in the low mass range, as can 
be seen by using the V-I vs J-K diagram shown in Fig. [9] 
where we present the color-color diagram of all objects; disk- 
less candidate members are used to estimate the cluster redden- 
ing by considering th e locus of expected photospheric colors by 
Kenv on & Hartmanr] d 1995b : we take advantage of the fact that 
in this diagram the expected photometric colors for low mass 
stars are almost vertical. 

In order to estimate the cluster reddening, we considered also 
the distribution of the J-K colors of diskless candidate members 
within fixed ranges of 0.2 magnitudes in the interval 1.3 <V- 
I< 3.0, and we considered, for each V-I bin, the peaks of the 
J-K distribution to derive the fiducial cluster sequence, indicated 
by the error bars in the bottom panel of Fig. [9] Comparison of 
this fiducial sequence with the lines allows us to confirm that the 
NGC 1893 mean reddening is E(B-V)=0.6±0.1. However, since 
the population of low mass stars is larger and spatially spread 
over the whole region affected by a larger differential reddening, 
we note that low mass candidate cluster members can individu- 
ally have E(B-V) larger or smaller values. As it will be shown 
later, by using IR colors we estimate that the reddening in this 
region is not larger than about E(B-V)=1.3 (Ay=4 mag). 

The spread of low mass candidate cluster members in the V 
vs. V-I diagram is not due only to differential reddening. In fact, 
if we select the sample of low mass diskless candidate members 
from the V-I vs J-K diagram with 0.5<E(B-V)< 0.7, we find 
that in the V vs. V-I diagram they show the same spread as all 
diskless and therefore we conclude that the low mass star spread 
in the V vs. V-I diagram could be also due to an intrinsic age 
spread and/or binarity. 

We note that the E(B-V) value we derived both for low 
and high mass candidate cluster members, is independent of 
other adopted parameters such as distance, ages and metallicity. 
However further spectroscopic observations should be obtained 
to derive spectral types and therefore individual reddening val- 
ues. 

6.3. Distance 

The distance of NGC 1 893 is a very uncertain parameter ranging 
in the literature from 3.2 to 6.0kpc (see Sect. [7). Using broad 
band optical or NIR photometry of candidate cluster members 
and the isochrone fitting method, it is very difficult to derive the 
cluster distance. In fact, very young clusters such as NGC 1893 
include a population of PMS stars, that lies in a ample region of 
the CMD diagram, and a small population, more or less rich, 
depending on the cluster age, of higher mass objects already 
on the MS phase and therefore on a well defined locus of the 
CMD. Nevertheless, the MS star distribution on the CMD is usu- 
ally almost vertical and therefore quite ineffective for an accu- 
rate distance derivation. In addition, for clusters younger than 
about 15 Myr, the flattening of the isochrones between the two 
phases (PMS and MS) ends on the MS at a given magnitude that 
is strongly dependent on the isochrone age. Unfortunately, un- 
less individual ages of PMS cluster members are well known, 
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Fig. 8. V vs. V-I diagram of all the objects (dots) with errors in V- 
I smaller than 0. 1 . Filled circles are the diskless candidate cluster 
members while squares indicate those with E(B-V) between 0.5 
and 0.7, a ccording to Fig.|71 Solid line is the solar metallicity 
ZAMS of Siess et all d2000t) at a distance of 3.6 kp c and redden- 
ing E( B-V)=0.4. Dashed lines indicate the 1.5 Myr Marigo et al. 
(2008) isochrone of solar metallicity for masses larger than 2 M Q 
at 3.6 kpc and for E(B -V)=0.5, 0.6 and . 7. The reddening vec- 
tor obtained using the iMunari & Carrarol JT9 96) reddening laws 
is also shown. 



it is very hard to constrain the cluster distance by looking for 
this magnitude, even if the cluster mean reddening is fixed and 
a large population of high and low mass cluster members is 
known. 

However, in the V vs. U-B diagram, the MS population is 
quite well separated from the bulk of the remaining objects that 
are field stars and candidate PMS cluster members. This is evi- 
dent in the upp er panels of F ig.fTOl obtai ned by using the UBV 
photometry by IMassev et all (1 19951) and ISharma et all (12007b . 
and where diskless candidate cluster members selected by us 
are indicated with filled circles. Among these objects, we indi- 
cate with squares those that according to the U-B vs. B-V dia- 
gram have 0.5<E(B-V)<0.7. We use these objects here to con- 
strain the cluster distance b y comparing again the solar metallic- 
ity isochrone of 1.5 Myr of Mari go et al.l (12008) reddened with 
E(B-V)=0.6. By taking advantage of the slightly curved shape 
of the isochrone in this CMD, we find that a good fitting is 
achieved if a distance of 3.6 kpc is adopted with an error of about 
200 pc. With our data, we can rule out for NGC 1893 a distance 
of 3200 p c or 4300 pc that are values previously given in the 
litera ture ( Sha rma et all 120071: iTapia et all 1 1 99 it IMassev et all 
H99l . 

A distance of 3.6 kpc for NGC 1893 is also consistent with 
that derived by using an independent method based on the 
nebula properties. In fact, as already done in previous works 
dPrisinzano et all l2005t iGuarcello et all 120071) . we assume that 
the nebula surrounding the cluster obscures many background 
objects and allows us to see mainly foreground stars, that are 
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Fig. 9. Color-color diagram of objects with error in color smaller 
than 0.1 (dots); solid line s are the photospheric colors by 
Kenv on & Hartmannl (1 19951) at the three E(B-V) values indi- 
cated on each line. Filled circles in the upper panel are the 
NGC 1893 candidate cluster members without circumstellar disk 
while error bars in the bottom panel indicate the cluster se- 
quence fiducial line derive as describ ed on the tex t . The red- 
dening vector obta i ned u sing the [Munari & Carrar ol (fT996h and 
Riek e & Lebofskvl (Il985l) reddening laws is also shown. 



mostly in MS, with distance smaller or equal to the nebula dis- 
tance. Among these latter, those located at the nebula distance 
define the blue envelope of the CMD diagram that can be fitted 
by the ZAMS locus. In the plausible assumption that the nebula 
is located at the same distance of NGC 1893 and by considering 
only very high reddened regions around the cluster (to minimize 
the spread due to background objects), it is possible to derive the 
cluster distance by fitting the ZAMS to the blue envelope of the 
CMD. 



By a visual inspection of our observations in the optical, 
JHK and Spitzer/IRAC bands, we note that the nebula is not 
evenly spatially distributed, lying mostly in the western region 
of our FoV. For this reason, we selected the MS stars to be 
used for the ZAMS isochrone fitting by using the most obscured 
regions. The latter have been selected by means of a redden- 



ngj 

in iDamiani et alj (120061 see this paper for details) by using H- 
K=0. 1 as limit for the intrinsic co lor expected for background 
giants. As in Damia nTet al.l (120061) . for the map computation we 
excluded all objects classified as candidate cluster members. The 
Ay map has been derived by using a cell size of 3.1'x 2.6' within 
the FoV defined by our NICS/JHK observations and by the E(K- 
H) map, obtained as the median of H-K colors in each cell minus 
the expected 0. 1 value. We find that the whole region is affected 
by 2.2 < Ay < 4, with the reddest regions being the western 
ones. As already mentioned, the whole region is not affected by 
very high Ay values and this can be also confirmed by the distri- 
bution of field stars in the J-H vs. H-K diagrams. This is likely 
due to a low absolute reddening affecting the galactic anti-center 
direction. 

By using the resulting Ay map, we selected all the objects 
within the regions with Ay > 3.8 and error in the V-I colors 
smaller than 0.1; the CMD diagram of these objects (indicated 
with dots) is shown in the bottom panels of Fig.[lO] where filled 
circles are diskless candidate members. The dashe d curve in 
the fig ure represents the solar metallicity ZAMS of Sies s et al.l 
(2000) shifted at the distance d=3.6kpc and using E(B-V)=0.4, 
that is the minimum reddening we associate to field stars ac- 
cording to the bright stars in this diagram, while the solid fine 
is the solar metallicity 1.5 Myr isochrone of Marigo et al. (2008) 
for stars with mass M>1.5M , adopting the same distance and 
E(B-V)=0.6, that is the mean cluster reddening we derived in the 
previous section. 

As already mentioned, we derive the cluster distance as that 
value for which a good ZAMS fitting to MS stars of the blue en- 
velope of this diagram is achieved. We note that, since the max- 
imum Ay is about 4, even in the most obscured regions, a little 
fraction of faint background objects could be visible. In fact, as 
shown in Fig [10] (bottom panel), objects with V> 21 do not fol- 
low the ZAMS shape, likely because this region is populated by 
field stars or extragalactic objects, that are expected to be visible 
in the galactic anti-center direction where the cluster NGC 1893 
is located. Therefore, for the ZAMS fitting we consider only field 
stars with V< 21; nevertheless, since the bright population of 
field stars (V< 17) is poorly populated, we considered, in this 
range, the candidate cluster members with V < 16 that are on 
the MS phase and with 0.5<E(B-V)<0.7, indicated with squares. 
We note that 16 is the magnitude of the Turn On, i.e. the 
magnitude where the PMS joins the cluster MS. We fitted the 
candidate MS cluster members assuming the mean cluster red- 
dening E(B-V)=0.6 and we found that a good fit is obtained both 
for the MS stars in the blue envelope and for the bright candi- 
date cluster members on MS if a distance of ~3.6kpc (see bot- 
tom panel of FigJTOb is used. By taking the solar galactocentric 
distance of 8.5 Kpc, we deduce that NGC 1893 is located at a 
distance of about 12.1 Kpc from the Galactic Center. 

7. Comparison with literature cluster parameters 

The cluster distance we derived (3.6±0. 2kpc) is marginally 



consistent with that recently derived by Sharma et al. (2007 



(d=32 50pc) by fitting the 4 Myr isochrone of Bertelli et al.l 
(119941) to the MS field stars and assuming E(B-V)=0.4, that is 
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Fig. 10. V vs . U-B diagram from the iMassev etaO d 1995b and 
Sharmaet al. (2007) photometry (upper panels) and V vs. V-I 
diagram of all the objects with errors in V-I smaller than 0. 1 and 
located in the regions with Ay > 3.8 (bottom panels). Filled 
circles are the diskless candidate cluster members; squares are 
those with 0.5<E(B-V)<0.7; solid line is the 3.0 Myr isochrone 
of solar metallicity using E(B-V)=0.6 while dashed line is the 
solar metallicity ZAMS of lSiess et all d2000h using E(B-V)=0.4; 
both are drawn assuming a distance of 3.6kpc. The reddening 
vector obtained using the iMunari & Carrar oT (fl996h reddening 
laws is also shown. 



the value taken for the MS field stars. This method is similar 
to that used b y us, since we note that the 4 Myr isochrone of 
iBertelli et all (11994 corresponds rather to the ZAMS locus, if 
we consider that stars 4 Myr old with masses smaller than about 

2 M Q are yet in PMS . 

We find a good agreement with the value derived by|Cuffev 
(fl97l who, using UBV magnitudes, derived that the redden- 
ing in front of the cluster amount to E(B-V)=0.4, while the 



distance is 3.6 kpc. Among the fi r st dist ance values published 
for NGC 1893. iBecker & Fenkartl (Il97lh found d=3700pc and 
Ay = 1.68, corresponding to E(B-V)= 0.54, by usin g simulta- 
neously two color-magnitude diagrams. Moffat (1972) derived 
a distance of 3980 pc and a reddening E(B-V)=0.55 using pho- 
tographic UB V photometry and the ZAMS fitting to the unred- 
dened CMD. Hump hrevsl (Il978l) found a distance modulus of 
12.52 corresponding to a distance of about 3200pc, by using 
V, B-V photometry of three O type main sequence stars of 
NGC 1893 in a paper devoted to study fundamental properties 

of the most luminous stars in our Galaxy. 

A larger distance of 4300 pc has been derived by Tapia et al] 
(Il99ll) by using Stromgren ubvy, Wfi and JHK photometry. They 
derived spectral types and approximate luminosity classes from 
the m\ and c\ indices that together with literature slit spectral 
types are used to estimate reddening and absorption and there- 
fore intrinsic colors and magnitudes for candidate cluster mem- 
bers. The corrected distance modulus is derived as the difference 
between the unabsorbed V magnitude and the absolute magni- 
tude predicted for the given spec tral type. A similar v alue, equal 
to 4800 pc, has been found by iFitzsimmonsI (JT993), by using 
again Stromgren ubvy photometry and the theoretical ZAMS fit- 
ting. 

CCD U BV photomet r y and multiobject fiber spectroscopy 
are used by Massev ~et al.1 (1 19951) to determine distance and red- 
dening equal to 4400pc and E(B-V)=0.53±0.2, respectively. 
These values are derived by using individual reddening from 
spectroscopy and appropriate spectral type-My calibration for 
O and B-type stars. The derived distance has been obtained by 
including only stars with inferred color excess E(B-V) between 
0.4 1 and 0.70 

IVallenari etal] d 19991) use J and K photometry and UBV lit- 
erature photometry to derive E(J-K) from 0.15 to 0.5, corre- 
sponding to 0.25<E(B-V)<0.8 in the whole region they study; in 
addition they find two main dark regions where E(J-K) is 0.45- 
0.55 (or E(B-V)=0.7-0.9) that is higher than the value derived in 
the regions outside the dark clumps. 

The most discordant distance value for NGC 18 93 with re- 
spect to that found in this paper is that derived by iMarco et al.l 
(1200 ll) equal to 6000 pc; they use photometric Stromgren in- 
dices (magnitude limit V^15.9) to compute the interstellar red- 
dening E(b - y) of cluster members. They find values up to 
E(b - y) =0.607 but they estimate the average interstellar red- 
dening E(b - y) = 0.33 ± 0.03 by considering only stars with 
E(b - y) < 0.4. They use this average reddening to derive indi- 
vidual magnitudes and colors ( Vo, (b - y)o and co). By fitting an 
empirical ZAMS in the My vs. en diagram they derive a dered- 
dened distance modulus of 13.9. Considering that the cluster 
is affected by a significant differential reddening, we suppose 
that th e photometric intrinsic parameters derived bv lMarco et al.l 
(1200 lh . based on the average reddening could be the origin of 
their overestimate of the cluster distance. 

8. Masses and Ages 

Assuming the cluster parameters derived by us in the previous 
sections, we was able to compute masses and ages of all ob- 
jects we selected as candidate cluster members. However, as we 
mentioned before, the position in the color-magnitude diagram 
of class II candidate members can be influenced by effects due 
to the presence of circumstellar disk. Since our data do not allow 
to disentangle the purely photospheric spectrum of the stars, to 
which theoretical models are referred, from that in excess due to 
such effects, we assume that masses and ages of diskless stars 
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Fig. 12. Age distribution of diskless and class II candidate mem- 
bers in the PMS region. Class II YSOs with anomalous blue col- 
ors are not included in this sample. 



can be accurately derived within the uncertainties due to pho- 
tometric errors, reddening and to the theoretical models, while 
masses and ages of class II stars should be taken with caution, 
especially those of the 170 candidate members falling outside 
the PMS region with ages apparently older than lOMyr. 

Masses and ages are co mputed by interpolating the solar 
metallicitv lSiess et al.l (12000) tracks on the V vs. V-I plane, us- 
ing the TRIGRI D idl function. F igurefTTl shows the V vs. V-I 
diagram with the lSiess et al. (2000) tracks and isochrones super- 
posed. 

For the diskless sample we find stars with masses between 
0.3 and 6.7 MgU we have further 1 1 objects classified as disk- 
less candidate members with mass larger than 7 M ; the median 
value and the standard deviation of diskless member ages are 
1.4Myr and 1.8 Myr, respectively. For the class II sample, if we 
discard the 170 stars with optical anomalous photometry, we find 
stars with masses between 0.2 and 6.9 M ; we have further 10 
objects classified as class II candidate members with mass larger 
than 7 M ; the median value and the standard deviation of ages 
of candidate members with disk are 1 .6 Myr and 2.2 Myr, respec- 
tively, with again few stars being 9- 10 Myr old. We conclude 
that, apart from the 170 objects with anomalous optical indices, 
the remaining class II stars show mass and age distributions very 
similar to those of diskless YSOs as it is shown in Fig.fTSIwhere 
the age distribution of diskless and class II candidate members 
are shown. This reflects the similar spread observed in the V 
vs. V-I diagram by the two populations of candidate members. 
Apparently the colors of class II candidate members are not sub- 
stantially modified by phenomena related to accretion. The peak 
in the observed mass function is among 0.4 and 0.5 M ; how- 
ever the IMF will be computed in an other work where we will 
discuss the biases due to the adopted methods and where we will 
compare the IMF with that of other clusters of similar age. 



3 the upper limit is imposed by the adopted ISiess et 23 HOOO) mod- 
els. 



9. Summary and conclusions 

We used new deep optical NIR data in the VRIJHK and H a 
bands together with literature X-ray and Spitzer-IRAC data to 
compile a multiband photometric catalog in the direction of the 
young distant open cluster NGC 1893, located towards the galac- 
tic anti-center. Using these data and different membership crite- 
ria, we made the most complete census of candidate members 
available in the literature for this cluster. Reddening-independent 
indices, involving optical and/or NIR colors, allowed us to de- 
rive a list of 1061 candidate members with a circumstellar disk, 
with 125 among them being also H„ emitters; among the disk 
bearing candidate members, 170 show anomalous optical mag- 
nitudes and colors that are likely due to effects of gas and/or dust 
in the disk. In addition, X-ray detections and optical magnitudes, 
allowed us to distinguish 415 candidate members without disk. 

We used these last 415 YSOs and their photometric prop- 
erties, to assess the cluster parameters, viz. interstellar redden- 
ing and distance. Disk-less candidate members are, in fact, the 
most reliable objects to be compared with theoretical tracks and 
isochrones, since their colors are purely photospheric and do not 
suffer from additional effects due to the presence of circumstel- 
lar disks. 

Bright diskless members in the U-B vs. B-V diagram, ob- 
tained from literature data, are consistent with an interstellar red- 
dening E(B-V)=0.6±0.1. This value is also consistent using the 
same bright candidate members in the independent V vs. V-I di- 
agram and most of the low mass diskless candidate members in 
the V-I vs. J-K diagram. However, we find evidence of differen- 
tial reddening in this region. 

Using the V vs. B-V (from the literature) diagram, and as- 
suming the mean cluster reddening E(B-V)=0.6, we find that 
the main sequence members selected in this work are well fitted 
by a 1.5 Myr isochrone, if a distance of 3.6±0.2kpc is adopted. 
We note that in the MS cluster region, isochrones degenerate 
and therefore our distance estimate is independent of the age 
adopted for MS members. The value of 3.6 kpc is consistent with 
the cluster distance we derive by using the main sequence blue 
envelope traced in the V vs. V-I diagram by foreground field 
stars located at the cluster distance. Therefore, using indepen- 
dent methods and data, we are able to firmly evaluate the cluster 
distance for which values between 3250 and 6000 pc are given 
in the most recent literature. 

Finally, we derived masses and ages of selected candidate 
members in the optical PMS region down to about 0.2 M ; me- 
dian ages for diskless and class II stars, are, respectively, 1 .4 and 
1 .6 Myr, with a standard deviation of 1 .8 and 2.2 Myr for the two 
samples. Few cluster candidate members are older with ages up 
to about 10 Myr. 

If we consider that 1001 of the 1061 class II candidate mem- 
bers are within the Chandra FoV, where we have selected the 
diskless candidate members, we find a disk fraction of about 
7 1 % of stars. This value is an upper limit since the diskless mem- 
ber selection is based on the X-ray observations that reach a mass 
limit higher than that we can reach with NIR observations used 
to select class II stars. However, the disk fraction we find is sim- 
ilar to the 67% value we found in Paper I and in a greement with 
the d isk fraction found in clusters of similar age (lHaisch et all 
2001), if we assume for NGC 1893 a median age of 1.5Myr. 
A detailed analysis of the disk fraction as a function of stellar 
masses and cluster location will be the subject of a forthcoming 
paper (Sanz-Forcada et al. 2010, in preparation). 

We therefore conclude that, despite its peculiar location in 
the Galaxy, NGC 1893 includes a rich population of young stars, 
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Fig. 11. V vs. V-I diagram of all objects with optical magnitudes in our cat alog (dots), diskless and class II candidate members are 
indicated with filled circles and empty squares, respectively. Lines are the lSiess et alJ (120001) tracks and isochrones of masses and 
ages indicated on each line. 



with general properties similar to those found in young clusters 
in the solar neighborhood. 

Acknowledgements. Part of this work was financially supported by the 
PRIN-INAF (P.I. Lanza) and the EC MC RTN CONSTELLATION 
(MRTNCT2006035890). We thank Donata Randazzo for carefull reading of the 
paper. 

References 

Baffa, C, Comoretto, G., Gennari, S., et al. 2001, A&A, 378, 722 
Becker, W. & Fenkart, R. 1971, A&AS, 4, 241 



Bertelli, G., Bressan, A., Chiosi, C, Fagotto, F, & Nasi, E. 1994, A&AS, 106, 
275 

Caramazza, M., Micela, G., Prisinzano, L., et al. 2008, A&A, 488, 211 
Cuffey, J. 1973, AJ, 78, 747 

Cutri, R. M., Skrutskie, M. E, van Dyk, S., et al. 2003, 2MASS All Sky Catalog 
of point sources, ed. U. of Massachusets and Infrared Processing and Analysis 
Center, IPAC/California Institute of Technology 

Daflon, S . & Cunha, K. 2004, ApJ, 617,1115 

Damiani, E, Prisinzano, L., Micela, G., & Sciortino, S. 2006, A&A, 459, 477 
Fitzsimmons, A. 1993, A&AS, 99, 15 

Flaherty, K. M., Pipher, J. L., Megeath, S. T., et al. 2007, ApJ, 663, 1069 
Gaze, V. F. & Shajn, G. A. 1952, Izvestiya Ordena Trudovogo Krasnogo 
Znameni Krymskoj Astrofizicheskoj Observatorii, 9, 52 



16 



L. Prisinzano et al.: Star formation in the outer Galaxy: the young open cluster NGC1893 



Ghinassi, F., Licandro, J., Oliva, E., et al. 2002, A&A, 386, 1 157 
Guarcello, M. G., Damiani, F., Micela, G, et al. 2010, ArXiv e-prints 
Guarcello, M. G, Micela, G, Damiani, F., et al. 2009, A&A, 496, 453 
Guarcello, M. G, Prisinzano, L., Micela, G, et al. 2007, A&A, 462, 245 
Gullbring, E., Hartmann, L., Briceno, C, & Calvet, N. 1998, ApJ, 492, 323 
Haisch, Jr., K. E., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 
Hillenbrand, L. A. 1997, AJ, 113, 1733 

Hiltner, W. A. 1966, in IAU Symposium, Vol. 24, Spectral Classification and 
Multicolour Photometry, ed. K. Loden, L. O. Loden, & U. Sinnerstad, 373-+ 

Hoag, A. A., Johnson, H. L., Iriarte, B., et al. 1961, Publications of the 
U.S. Naval Observatory Second Series, 17, 345 

Humphreys, R. M. 1978, ApJS, 38, 309 

Kenyon, S. J. & Hartmann, L. 1995, ApJS, 101, 117 

Landolt, A. U. 1992, AJ, 104, 340 

Leggett, S. K., Currie, M. J., Varricatt, W. P., et al. 2006, MNRAS, 373, 781 

Marco, A., Bernabeu, G, & Negueruela, I. 2001, AJ, 121, 2075 

Marco, A. & Negueruela, I. 2002, A&A, 393, 195 

Marigo, P., Girardi, L., Bressan, A., et al. 2008, A&A, 482, 883 

Massey, P., Johnson, K. E., & Degioia-Eastwood, K. 1995, ApJ, 454, 151 

Moffat, A. F. J. 1972, A&AS, 7, 355 

Munari, U. & Carraro, G. 1996, A&A, 314, 108 

Prisinzano, L., Damiani, F., Micela, G, & Sciortino, S. 2005, A&A, 430, 941 
Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 

Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 

2006, ApJS, 167, 256 
Rolleston, W. R. J., Brown, P. J. E, Dufton, P. L., & Fitzsimmons, A. 1993, A&A, 

270, 107 

Rolleston, W. R. J., Smartt, S. J., Dufton, P. L., & Ryans, R. S. I. 2000, A&A, 
363, 537 

Sanchez, S. F, Aceituno, J., Thiele, U., Perez-Ramirez, D., & Alves, J. 2007, 
PASP, 119, 1186 

Sharma, S., Pandey, A. K., Ojha, D. K., et al. 2007, MNRAS, 380, 1141 
Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 
Stetson, P. B. 1987a, PASP, 99, 191 
Stetson, P. B. 1987b, PASP, 99, 191 
Stetson, P. B. 1990, PASP, 102, 932 
Stetson, P. B. 1994, PASP, 106, 250 
Stetson, P. B. 2000, PASP, 112, 925 
Stetson, P. B. 2005, PASP, 117, 563 

Tapia, M., Costero, R., Echevarria, J., & Roth, M. 1991, MNRAS, 253, 649 
Vallenari, A., Richichi, A., Carraro, G, & Girardi, L. 1999, A&A, 349, 825 
Wilson, T. L. & Matteucci, F. 1992, A&A Rev., 4, 1 

Wouterloot, J. G. A., Brand, J., Burton, W. B., & Kwee, K. K. 1990, A&A, 230, 
21 

Zinnecker, H. & Yorke, H. W. 2007, ARA&A, 45, 481 



L. Prisinzano et al.: Star formation in the outer Galaxy: the young open cluster NGC1893, Online Material p 1 



Appendix A: Photometric calibration 

A.1. DOLORES 

Photometric calibration of DOLORES observations was com- 
puted by using the images of the Landolt (1992) standard field 
SA98. We used the v, r, i, h a instrumental magnitudes, and con- 
sidered the V, R and I magnitudes of the Johnson-Kron -Cousins 
photo metric system of standard stars in these fields dStetsonL 
2000); H ff is the magnitude in the corresponding non-standard 
filter. We calculated the transformation coefficients between the 
instrumental an d standard systems by using the CCDSTD code 
(Stetson, 2005) and the following equations: 

v = v + A +A l Q+A 2 {V-I)+A 3 (V-I) 2 +A 4 XY + 

+A 5 X + A 6 Y + A-, X 2 + A 8 Y 2 
r = R + B + BiQ + B 2 {R- I) + B 3 (R- I) 2 + B 4 XY + 

+B 5 X + B 6 Y + B 7 X 2 + B S Y 2 (A.l) 
i = i + q ) + C 1 Q + C 2 (V-I) + C 3 (V-I) 2 + C 4 XY + 

+C 5 X + C 6 Y + d X 2 + C 8 Y 2 
K = H a + D Q + D l Q + D 2 (R-I) + £> 3 (R- I) 2 + D 4 XY + 

+D 5 X + D 6 Y + D 7 X 2 + D 8 Y 2 

where Q is the airmass, Aq, Bq and Co are the magnitude zero 
points; A\, B\ and C\ are the extinction coefficients; A 2 , B 2 and 
C 2 and A3, B3 and C3 are the color terms, and the rest of terms are 
related to the geometrical coordinates (X, Y) of the plate. Stars 
with magnitude residual values above the 3cr level were not con- 
sidered for the coefficient fitting. Coefficients calculated during 
calibration are listed in Table lA. ll DOLORES extinction terms 
are set to those provided by the telescope web site. We followed 
the same nomenclature for the Hff filter, but since we have no 
standard measurements in this band, we applied the same cal- 
ibration as for the R band (therefore D, = C,) in the case of 
DOLORES. 

A.2. CAFOS 

We tried to calibrate the CAFOS observations following the 
same procedure as for DOLORES, using the site extinction 
terms as provided by Sanc hez et alj d2007l) . However we no- 
ticed substantial discrepancies in the final magnitudes calculated 
for CAFOS when compared to DOLORES results. The vari- 
able atmospheric conditions detected during the observations, 
also noted in the seeing, could be responsible of the inconsis- 
tencies. We decided to use the DOLORES photometry to cal- 
ibrate the CAFOS fields: we selected the local standards used 
in DOLORES fields that are also common to the CAFOS FoV, 
and used a total of 417 CAFOS counterparts as standards for 
the calibration, so we calibrated the CAFOS observations in the 
same photometric system of DOLORES. The second and third 
nights of CAFOS observations also required to set the "cloud" 
variable in the CCDSTD setting. As explained bv lStetsonl d2000t) 
this variable allows to consider different atmospheric conditions 
(thin clouds) during the night, setting all the extinction terms 
to in this case (Table lATT i. and fitting the zero points and the 
other terms. CAFOS observations used the DOLORES standard 
stars falling in the field to calibrate the Ha observations, with D, 
coefficients as listed in Table lA. ll 



images obtained with SNAP for each filter (JHK) and for each of 
the 16 fields observed around the cluster NGC 1893. For most of 
the images we performed the PSF photometry using the Moffat 
function with B = 2.5, typically used to model stellar profiles 
as an analytic first approximation to the PSF, while for 6 im- 
ages we used the more complex Penny function in order to take 
into account the elongated shape of the PSF likely due to some 
aberration. In addition we treated the variable PSF by using the 
DAOPHOT option that considers a PSF which varies quadrati- 
cally with position in the frame. 

Since the PSF photometry is relative to the model stellar pro- 
file of a given frame, we needed to derive the aperture correction 
to the instrumental photometry; to this aim we performed aper- 
ture photometry at different radii on a sample of isolated and rel- 
atively bright stars fro m which we de rived growth curves with 
the DAOGROW code dStetsonL Il990h . The aperture correction 
has been computed using the selected stars as the median of the 
difference between the PSF magnitudes and the aperture magni- 
tudes obtained at the radius including all the stellar flux. 

In order to discard false identifications due to the spots ei- 
ther of very bright objects or just falling at the edge of the 
images, we selected the list of objects detected by consider- 
ing only those following the typical exponential profile of the 
magnitude errors. Then for each of the 16 observed fields , we 
merged with DAOMATCH/DAOMASTER dStetsonL fl987bl) the 
three lists with JHK magnitudes in order to have a single list with 
objects having at least two among the JHK magnitudes. 

The NICS filters used for our observations (Js, H, K') are 
those o f the Mauna Kea Obs ervatories (MKO) near-infrared fil- 
ter set (Ghinassi et al., 2002) and therefore we calibrate our cat- 
alog in the MKO photometric system. To this aim we use as 
standard stars, the 2MASS counterparts falling in our fields with 
PH.QUAL flag equal to 'AAA' and CC_FLG flag equal to '000' . 
We first converted the 2MASS catalog in the MKO system by us- 
ing the inverted transformations giv en in the 2MASS w eb pagf] 
maintened by J. Carpenter (see also lLeggett et al.Ll2006l) . 

By using instrumental magnitudes for the standard stars, that 
we indicate as j, h, k and those of the 2MASS catalog in the 
MKO system, that we indicate as Jmko, Hmko, Kmko we per- 
formed for each field a linear fit of the magnitude differences as 
a function of the colors as in the following equations: 



(j - Jmko) = A + Ai(j -h) 
(j-JMKo) = B + B 1 (j-k) 
(h-HMKO^Co + QCj-h) (A.2) 
(h-HMKO^Do+DKh-k) 
(k — Kmko) = E +Ei(j -k) 
(k-K M Ko) = F +F 1 (h-k). 

From the linear fit we derived the zero points, indicated by coef- 
ficients with subscript and the color terms, indicated by coef- 
ficients with subscript 1 given in Table lA. 2l The final photomet- 
ric catalog in the MKO system has been obtained by imposing 
these coefficients and the analogous equations to the instrumen- 
tal magnitudes of all detected objects. 



A.3. NICS 



We use d the standalone DAOPHOT II/ALLSTAR code dStetsonL 
1987b) to obtain the instrumental photometry for the combined 



4 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4b.html 
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Table A.l. Coefficients of the transformations to the standard system, given in the Eq.s lA.ll for each filter and night for DOLORES 



(D)° and CAFOS (C) 


Instr. Date 


Filter 


Coeff. 




A u 


B„ C u D 


coefficient indices 

















1 


2 


3 


4 


5 


6 


7 


8 


D 2007/09/21 


V 


Aj 


-0.874 


0.15 


0.087 






-0.068 


0.138 




-0.045 








±0.014 




±0.007 






±0.006 


±0.021 




±0.010 


D 


R 


B, 


-0.823 


0.11 


-0.002 






-0.072 


-0.007 












±0.008 




±0.009 






±0.004 


±0.003 






D 


I 


Q 


-0.531 


0.07 


-0.035 






0.420 


0.393 


-0.191 


-0.181 








±0.020 




±0.006 






±0.026 


±0.017 


±0.010 


±0.008 


D 2007/10/18 


V 


At 


-1.039 


0.15 


0.103 






-0.007 


0.015 




-0.006 








±0.004 




±0.002 






±0.002 


±0.007 




±0.003 


D 


R 


B, 


-1.100 


0.11 


0.015 






0.013 


0.011 












±0.003 




±0.004 






±0.002 


±0.002 






D 


I 


Ci 


-0.572 


0.07 


0.018 






0.237 


0.224 


-0.081 


-0.117 








±0.009 




±0.004 






±0.019 


±0.010 


±0.011 


±0.005 


D 2007/11/14 


V 


Ai 


-1.019 


0.15 


0.090 






-0.014 


0.072 




-0.026 








±0.013 




±0.004 






±0.016 


±0.011 




±0.005 


D 


R 


B, 


-1.050 


0.11 


-0.007 






0.043 


0.002 


-0.026 










±0.005 




±0.004 






±0.008 


±0.002 


±0.003 




D 


I 


Q 


-0.537 


0.07 


-0.019 






0.307 


0.239 


-0.136 


-0.128 








±0.010 




±0.004 






±0.014 


±0.010 


±0.006 


±0.005 


C 2007/10/11 


V 


Ai 


2.482 


0.118 


-0.022 


0.015 


-0.026 


-0.236 


-0.127 


0.139 


0.086 








±0.019 


±0.004 


±0.015 


±0.005 


±0.011 


±0.022 


±0.023 


±0.010 


±0.012 


C 


R 


B, 


1.820 


0.045 


0.030 


-0.056 


-0.147 


0.261 


0.151 












±0.042 


±0.011 


±0.049 


±0.031 


±0.061 


±0.082 


±0.091 






C 


I 


Q 


2.467 


0.051 


-0.082 


-0.015 


-0.039 


0.206 


0.344 


-0.132 


-0.208 








±0.021 


±0.005 


±0.009 


±0.004 


±0.025 


±0.030 


±0.040 


±0.021 


±0.022 


c 


H„ 


Di 


-1.344 


1.088 


0.340 




0.063 


-0.333 


-0.384 


0.205 


0.238 








±0.027 


±0.323 


±0.036 




±0.014 


±0.028 


±0.029 


±0.012 


±0.016 


C 2008/01/05 


V 


Ai 


2.642 




0.063 


-0.012 


-0.018 


0.190 


0.118 


-0.094 


-0.098 








±0.012 




±0.011 


±0.004 


±0.008 


±0.014 


±0.013 


±0.007 


±0.007 


C 


R 


B, 


1.940 




0.038 


-0.037 


-0.010 


0.201 


0.129 


-0.099 


-0.131 








±0.010 




±0.014 


±0.008 


±0.007 


±0.012 


±0.011 


±0.005 


±0.006 


c 


1 


Ci 


4.497 




-0.108 


-0.006 


-0.222 


-0.200 


0.233 


0.174 


-0.043 








±0.021 




±0.011 


±0.003 


±0.015 


±0.027 


±0.024 


±0.013 


±0.013 


c 


H„ 


Di 


2.948 




-0.373 


0.374 


0.027 


-0.266 


-0.594 


0.162 


0.347 








±0.105 




±0.154 


±0.127 


±0.060 


±0.130 


±0.100 


±0.052 


±0.059 


C 2008/01/09 


H„ 


Di 


1.724 




0.045 


-0.013 


0.118 


0.322 


-0.266 


-0.176 


0.028 








±0.045 




±0.051 


±0.035 


±0.026 


±0.049 


±0.050 


±0.019 


±0.024 



Notes: "DOLORES observations in Ha have same coefficients as for the R filter 
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Table A.2. Zero points (coefficients with subscript 0) and color terms (coefficients with subscript 1) derived from the linear fit of 
equations ( IA.2I ) to transform instrumental NICS magnitudes to the MKO system. 



10 



11 



12 



13 



14 



15 



16 



Ap 
-22.505 
± 0.013 
-22.577 
± 0.026 
-22.388 
± 0.025 
-22.410 
± 0.016 
-22.504 
± 0.012 
-22.471 
± 0.014 
-22.063 
± 0.010 
-22.627 
± 0.025 
-22.336 
± 0.014 
-22.409 
±0.016 
-22.358 
± 0.015 
-22.348 
± 0.025 
-22.469 
± 0.032 
-22.401 
±0.019 
-22.538 
± 0.021 
-22.489 
± 0.026 



0.039 
±0.018 

0.064 
± 0.040 
-0.014 
± 0.036 
-0.107 
± 0.028 

0.044 
± 0.019 

0.053 
± 0.022 
-0.023 
± 0.024 

0.257 
± 0.035 
-0.023 
± 0.032 
-0.032 
± 0.024 

0.057 
± 0.020 
-0.091 
± 0.038 

0.123 
± 0.050 

0.029 
± 0.032 

0.215 
± 0.032 

0.067 
± 0.041 



Bp 
-22.487 
± 0.007 
-22.548 
± 0.010 
-22.397 
± 0.009 
-22.459 
± 0.006 
-22.485 
± 0.005 
-22.450 
± 0.006 
-22.073 
± 0.005 
-22.493 
±0.011 
-22.346 
± 0.006 
-22.428 
± 0.007 
-22.332 
± 0.006 
-22.389 
±0.010 
-22.419 
±0.012 
-22.386 
± 0.006 
-22.436 
± 0.008 
-22.459 
± 0.008 



Bi 

0.027 
± 0.013 

0.049 
± 0.028 

-0.004 
± 0.026 

-0.076 
± 0.022 

0.034 
± 0.015 

0.048 
± 0.017 

0.011 
± 0.018 

0.198 
± 0.037 

0.010 
± 0.026 

-0.004 
± 0.017 

0.042 
± 0.013 

-0.079 
± 0.029 

0.110 
± 0.041 

0.011 
± 0.025 

0.176 
± 0.024 

0.058 
± 0.029 



Co 
-22.522 
± 0.015 
-22.740 
± 0.029 
-22.581 
± 0.029 
-22.519 
± 0.020 
-22.609 
±0.013 
-22.581 
± 0.015 
-22.074 
± 0.013 
-22.670 
± 0.029 
-22.450 
±0.017 
-22.557 
± 0.019 
-22.553 
± 0.018 
-22.536 
± 0.028 
-22.531 
± 0.039 
-22.551 
± 0.022 
-22.671 
± 0.022 
-22.554 
±0.031 



-0.056 
± 0.019 

0.014 
± 0.044 

-0.039 
± 0.041 

-0.231 
± 0.032 

-0.101 
± 0.020 

-0.053 
± 0.025 

-0.069 
± 0.029 

-0.067 
± 0.037 

-0.026 
± 0.037 

-0.070 
± 0.027 

0.005 
± 0.023 

-0.144 
± 0.041 

-0.072 
± 0.059 

-0.044 
± 0.035 

0.077 
± 0.033 

-0.113 
± 0.050 



Do 
-22.598 
± 0.015 
-22.673 
± 0.042 
-22.617 
± 0.048 
-22.757 
± 0.034 
-22.708 
± 0.023 
-22.583 
± 0.023 
-22.076 
± 0.015 
-22.664 
±0.051 
-22.260 
± 0.037 
-22.562 
± 0.020 
-22.507 
±0.012 
-22.646 
± 0.047 
-22.412 
± 0.063 
-22.661 
± 0.054 
-22.502 
± 0.032 
-22.565 
± 0.044 



Pi 

-0.108 
± 0.043 

0.142 
± 0.099 

-0.020 
± 0.098 

-0.255 
± 0.082 

-0.100 
± 0.058 

0.078 
± 0.060 

0.094 
± 0.052 

0.114 
±0.105 

0.580 
±0.106 

0.109 
± 0.050 

0.120 
± 0.031 

-0.034 
±0.108 

0.441 
±0.167 

-0.169 
±0.107 

0.272 
± 0.070 

0.138 
±0.106 



Ep 
-21.957 
± 0.009 
-22.178 
± 0.013 
-21.928 
± 0.010 
-22.079 
± 0.008 
-22.116 
± 0.007 
-22.081 
± 0.007 
-21.682 
± 0.006 
-22.145 
± 0.016 
-22.005 
± 0.007 
-22.019 
± 0.009 
-21.969 
± 0.008 
-22.049 
±0.013 
-22.031 
± 0.019 
-21.939 
± 0.010 
-22.038 
± 0.009 
-22.098 
±0.011 



-0.007 
± 0.014 

0.072 
± 0.034 
-0.010 
± 0.029 
-0.030 
± 0.025 

0.014 
± 0.017 

0.056 
± 0.017 
-0.004 
± 0.019 

0.246 
± 0.048 

0.074 
±0.031 

0.032 
± 0.020 

0.078 
± 0.013 
-0.012 
± 0.033 

0.044 
± 0.057 

0.080 
t 0.032 

0.077 
± 0.022 

0.070 
± 0.032 



Fo 
-21.964 
± 0.015 
-22.077 
± 0.047 
-21.942 
± 0.043 
-22.148 
± 0.035 
-22.128 
± 0.025 
-21.998 
± 0.021 
-21.683 
± 0.014 
-22.510 
± 0.064 
-21.918 
± 0.047 
-21.980 
±0.019 
-21.882 
±0.010 
-22.098 
± 0.052 
-21.936 
± 0.077 
-21.883 
± 0.064 
-21.966 
± 0.030 
-22.000 
± 0.041 



-0.015 
± 0.044 

0.205 
± 0.114 
-0.024 
± 0.089 
-0.157 
± 0.083 
-0.041 
± 0.063 

0.189 
t 0.056 
-0.001 
± 0.049 
-0.891 
±0.133 

0.242 
±0.132 

0.078 
t 0.049 

0.162 
± 0.028 
-0.110 
±0.121 

0.224 
± 0.207 

0.093 
±0.129 

0.133 
± 0.062 

0.203 
±0.101 



Table 6. Optical photometry of Class 0/1 and Class II candidate members 



Seq. 


XID 


RA(2000) 


Dec(2000) 


V 




R 




I 




Ha 


Num. 




[deg] 


[deg] 
















1 


55 


80.6220775 


33.5139642 


9.045± 


0.016 


9.241± 


0.047 






0.036 


2 


779 


80.7594680 


33.5270160 


11.132± 


0.014 


10.978± 


0.017 


10.817+ 


0.011 


0.009 


3 


664 


80.7418460 


33.4434210 


12.147+ 


0.014 


11.923+ 


0.050 


11.725+ 


0.016 


0.009 


4 


276 


80.6880800 


33.4065610 


12.269± 


0.010 


12.034± 


0.035 


11.793+ 


0.011 


0.014 


5 




80.7888715 


33.5006385 


12.692± 


0.022 


12.103+ 


0.007 


11.721+ 


0.013 


0.037 


6 




80.6358343 


33.5480553 


12.733± 


0.025 






1 1.819+ 


0.006 




7 


887 


80.7798330 


33.5475840 


12.938+ 


0.011 


12.491± 


0.008 


12.110+ 


0.014 


0.009 


8 


345 


80.6968720 


33.4187280 


13.047± 


0.011 






12.628± 


0.039 




9 


773 


80.7583743 


33.5198142 


13.084± 


0.015 


12.211+ 


0.016 


11.365+ 


0.021 


0.007 


10 




80.7442230 


33.5280269 


13.105+ 


0.007 


12.588+ 


0.025 


12.083± 


0.022 


0.006 



Notes - This table is published in its entirety in the electronic edition of Astronony & Astrophysics. A portion is shown here for guidance 
regarding its form and content. 
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Table 7. Optical photometry of diskless candidate members 



Seq. 


XID 


RA(2000) 


Dec(2000) 


V 


R 


I 




Num. 




[deg] 


[deg] 










1069 


246 


80.6833230 


33.4407100 


10.207± 0.042 




9.904+ 0.011 




1070 


281 


80.6885760 


33.3713500 


10.724± 0.009 


10.605+ 0.008 


10.254+ 0.020 


0.035 


1071 


266 


80.6865470 


33.4434930 


11.324+ 0.034 




10.978+ 0.013 




1072 


11 


80.5768170 


33.4725590 


11.716+ 0.019 




11.412+ 0.038 




1073 


5 


80.5541260 


33.5671010 


12.154+ 0.012 




11.487+ 0.025 




1074 


362 


80.6991250 


33.3685801 


12.451+ 0.006 


12.004+ 0.012 


11.398+ 0.014 


0.041 


1075 


305 


80.6923210 


33.4224650 


12.476+ 0.020 




11.913+ 0.022 




1076 


651 


80.7396750 


33.4333360 


12.710+ 0.016 


12.465+ 0.030 


12.180+ 0.019 


0.006 


1077 


420 


80.7068503 


33.4259343 


12.741+ 0.018 




12.196+ 0.011 




1078 


155 


80.6625410 


33.4408110 


13.123+ 0.022 




12.422+ 0.006 





Notes - This table is published in its entirety in the electronic edition of Astronony & Astrophysics. A portion is shown here for guidance 
regarding its form and content. 



Table 8. NIR photometry of Class 0/1 and Class II candidate members 



Seq. 


RA(2000) 


Dec(2000) 


J 


H 


K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


Num. 


[deg] 


[deg] 

















1 80.6220775 

2 80.7594680 

3 80.7418460 

4 80.6880800 

5 80.7888715 

6 80.6358343 

7 80.7798330 

8 80.6968720 

9 80.7583743 

10 80.7442230 



33.5139642 
33.5270160 
33.4434210 
33.4065610 
33.5006385 
33.5480553 
33.5475840 
33.4187280 
33.5198142 
33.5280269 



10.468+ 0.016 
11.781+0.043 
11.574+ 0.016 
10.905+ 0.039 
11.859+ 0.070 
11.403+ 0.036 
12.439+ 0.016 



12.814+ 0.089 
11.519+ 0.019 
10.434+ 0.029 
11.430+ 0.059 
11.630+ 0.048 
12.501+ 0.034 



10.583+ 0.027 
11.746+ 0.019 
11.306+ 0.008 
10.116+ 0.018 
11.477+ 0.061 
11.099+ 0.036 
12.326+ 0.013 



8.290+ 0.053 
10.326+ 0.041 
11.045+ 0.002 
11.331+ 0.002 
9.262+ 0.041 
10.686+ 0.001 
10.927+ 0.001 
11.992+ 0.016 
9.367+ 0.041 
10.697+ 0.002 



8.196+ 0.036 
10.219+ 0.001 
10.905+ 0.002 
11.255+ 0.002 
8.894+ 0.042 
10.747+ 0.001 
10.959+ 0.001 



9.376+ 
10.654= 



0.048 
. 0.002 



8.179+ 0.001 
10.099+ 0.004 
10.802+ 0.005 
11.373+ 0.017 
8.568+ 0.002 
10.726+ 0.005 
10.963+ 0.009 
11.971+0.036 
9.295+ 0.002 
10.517+ 0.017 



8.206+ 0.036 
9.718+ 0.015 
10.662+ 0.016 
11.399+ 0.051 
8.380+ 0.003 
10.811+ 0.025 
10.992+ 0.026 
11.845+ 0.079 
9.348+ 0.010 
10.015+ 0.054 



Notes - This table is published in its entirety in the electronic edition of Astronony & Astrophysics. A portion is shown here for guidance 
regarding its form and content. 



Table 9. NIR photometry of diskless candidate members 



Seq. 


RA(2000) Dec(2000) 


J 


H 


K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


Num. 


[deg] [deg] 

















1069 
1070 
1071 
1072 
1073 
1074 
1075 
1076 
1077 
1078 



80.6833230 
80.6885760 
80.6865470 
80.5768170 
80.5541260 
80.6991250 
80.6923210 
80.7396750 
80.7068503 
80.6625410 



33.4407100 
33.3713500 
33.4434930 
33.4725590 
33.5671010 
33.3685801 
33.4224650 
33.4333360 
33.4259343 
33.4408110 



10.041+ 0.021 
10.545+ 0.025 
11.348+ 0.035 
10.811+0.020 

11.485+ 0.017 
11.667+ 0.022 
12.128+ 0.094 
11.902+ 0.017 



9.394+ 0.053 
9.906+ 0.019 
10.508+ 0.018 
11.103+ 0.035 
10.577+ 0.027 

11.361+0.019 
11.442+ 0.023 
12.254+ 0.150 
11.711+0.018 



9.532+ 0.043 

10.628+ 0.012 
10.982+ 0.014 
10.545+ 0.010 

11.266+ 0.009 
11.349+ 0.012 
11.851+0.105 
11.628+ 0.010 



9.559+ 0.027 
9.840+ 0.039 
10.582+ 0.002 
10.957+ 0.001 
10.506+ 0.035 
10.259+ 0.070 
11.316+ 0.014 
11.177+ 0.002 
11.481+ 0.004 
11.494+ 0.001 



9.605+ 0.040 
9.914+ 0.001 
10.574+ 0.002 
10.946+ 0.001 
10.449+ 0.001 
10.395+ 0.001 
11.289+ 0.071 
11.143+ 0.002 
11.512+ 0.005 
11.490+ 0.002 



9.580+ 0.003 
9.903+ 0.002 
10.637+ 0.009 
10.954+ 0.006 
10.503+ 0.005 
10.375+ 0.004 
11.259+ 0.030 
11.160+ 0.007 
11.302+ 0.049 
11.474+ 0.009 



9.658+ 0.010 
9.924+ 0.009 
10.697+ 0.029 
11.029+ 0.023 
10.539+ 0.017 
10.370+ 0.014 
11.136+ 0.062 
11.102+ 0.022 
11.431+ 0.040 
11.497+ 0.023 



Notes - This table is published in its entirety in the electronic edition of Astronony & Astrophysics. A portion is shown here for guidance 
regarding its form and content. 



